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SUMMARY 

The computational formulas used in the Brouwer theory of 
an artificial satellite are modified by a method similar to that 
suggested by R. H. Lyddane to remove the singularities at cir- 
cular and equatorial orbits. The Brouwer equations are also 
modified to reflect the use of the Vinti form for the force func- 
tion and to show all factors of e" and sin i" explicitly. 
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MODIFICATION OF BROUWER'S SOLUTION 
FOR ARTIFICIAL SATELLITES TO INCLUDE 
SMALL ECCENTRICITIES AND INCLINATIONS 

INTRODUCTION 

As Brouwer (Reference 1) states, the singularities in his formulas for 
small eccentricities and small inclinations are apparent since singularities do 
not exist in the coordinates. However, for numeric evaluation these singulari- 
ties are quite real and cause erroneous results when the eccentricity or inclin- 
ation is very small. Brouwer suggests that in these singular cases the formulas 
be modified to obtain expressions for the perturbations in coordinates. Although 
this approach is feasible the modifications would be quite extensive and give 
rise to different algorithms for different orbits. Lyddane (Reference 2) has 
suggested a modification which is not too far removed from Brouwer's equations 
and yields a single algorithm for all orbits. 

Although the primary purpose of the present modification is to remove 
singularities at small eccentricities and inclinations, Brouwer's formulas are 
also modified to meet the following requirements: (1) To show factors of e" and 
sin i " explicitly rather than implicitly so that limits become obvious as either 
of these values approach zero, (2) To change the form of the force function to 
the form used by Vinti as Brouwer recommends and (3) To adapt the formulas 
so that they are better suited for machine calculations. In view of this last 
requirement it is noted that many evaluations of quadratic polynomials are 
required in computing the perturbations. Hence, a single function for generating 
quadratic evaluations (such as a macro instruction) will simplify the coding of 
the procedure. For this reason all quadratic evaluations are shown separately. 

The development of the modified formulas is one of straightforward algebraic 
manipulation but tedious requiring extreme caution to avoid errors. The formulas 
contained herein have been verified several times algebraicly by independent 
checks and verified numerically by programming the procedure and comparing 
the numeric results with existing computer programs using the formulas of 
Brouwer. To avoid the possibility of typographical errors the procedure was 
programmed for a computer using the formulas in the review copy and the 
results of this program were compared with previous results. 


METHOD 

The basic modification to Brouwer's formulas is merely to replace the 
classical Keplerian elements a, e, i, h, g, l by simple functions of these elements 
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which are nonsingular at zero eccentricity or inclination. Many such functions 
exist, however, we have chosen the following which vary slightly from those used 
by Lyddane: a, A. - Z + g + h, /Zj = e cos Z , /a 2 = e sin Z , Vj - sin i cos h, and 
v 2 = sin i sin h. The osculating values of these functions may be obtained by 
Taylor series expansions about the mean values and ignoring second and higher 
order terms. Thus, 


a - a " + 8 a 


k - k" +8Z+8g+8h 
= (e" + Se) cos Z " - (e"8Z)sinZ" 

p. 2 = (e"+8e)sinZ" + (e" 8Z ) cos Z " 

Vj = cosh"(sini" + 8icosi") - (sini"8h)sinh" 
v 2 = sinh"(sini" + 8 i cos i " ) + ( s in i " 8h) cos h" 


where the operator 8 represents the sum of Brouwer's long and short period 
terms. As Lyddane points out, the higher terms of the Taylor series are singu- 
lar and ignoring them is mathematically unjustifiable. However, the Taylor 
series approach yields the same results as those obtained by Lyddane. Brouwer 
uses l ' and g' in the computation of the short period terms indicating that the 
Z" and g" might be used. In the present case l" and g" must be used since Z ' 
and g' may be ill-defined. 


FORCE FUNCTION AND BASIC CONSTANTS 
Using the form 


U 



11 Jk ( 7 ) P k (sin/3) 

k = 2 


2 



as the adopted force function and the basic constants 


a 0 ° a 


e Q - e 


x o - 1 


go 


h o + go + ^0 


R 


semi-major axis constant, 
eccentricity constant, 
inclination constant, 

right ascension of ascending node constant, 
argument of perigee constant, 
mean anomaly constant, 

1/2 -3/2 

M a o > 

mean longitude constant, 
equatorial radius, 


and comparision with Brouwer's form for the force function gives 


1 2 

k 2 = 2 JjR ! , 

k 4 = - |j 4 E 4 , 

The Brouwer abbreviations become 



*3.0 = -J 3 R 3 . 

*5.0 = - J S RS • 

8 - cos i Q , 


3 



r 2 

= 5 J = (1) ’ 

y 3 

~ Ja 


y 4 

■ 

y s 

= - Js 

(f.y 

yi 

= 2 J ! (f ) ■ 

yi 

= " Ja 

(ti 

yi 

3 /R\ 4 

8 ^4 J > 

yi 

= - Js 

(IT 


where 

P = a 0 V 2 ■ 

/ 

The y i and y i do not appear in the present development these being replaced by 
the Jj and the abbreviation y - - 1/2 (R/p). 

SECULAR TERMS 


PJI 

25 t } 2 + 16 77 

- 15 , 

P i 2 

- 90 17 2 ~ 96 77 + 

30 , 

P 13 

- 25 77 2 + 144 77 

+ 105 

q u 

= 3 6* 2 - 1 , 


q J 2 

= p I3 e 4 + p I2 

tf 2 + p 

11 > q l3 

= 3e 0 2 (35 (9 4 - 30 <9 2 

0 ) 


A l = Tj Jj 

2 q ll + 

%y 2 (J 2 2 q i 2 

“ SJ 4 q i3 )J . 


P«I 

- 25 77 2 + 24 77 

- 35 , 

Pg 2 

= - 126 17 2 - 192 7 J 

+ 90 


4 


Pg3 

= 45 rj 2 + 360 t] + 385 , 

q = se 2 - 1 , 

Og 2 

= P g3 + P g2 ° 2 + P g i • 

Pg 4 = - 9r , 2 + 21 , 

Pgs 

= 126 t] 2 - 270 , 

P g6 = - 189 r ? 2 + 385 , 


Og3 = P. 6* 4 

+ p - P . 

*g5 *g 4 


A g = J 2 0 + s y 2 

(j 2 2 O g 2 -5J 4 q g3 ) 

Phi 

- 9 77 2 + 1 2 77 — 5 , 

Ph 2 “ 5 77 2 + 3677 + 35 

o h i 

= - 2 , 

Oh 2 “ “ 4 (Ph 2 ^ 2 _ Phl) 


Oh 3 = 4 ( 37 ? 2 

- 5) (7 (9 2 — 3) , 


A h ~ ^ _J 2 Oh 1 + 8 

7 2 (J2 2 1 h! -5J 4 q,,3)] 

1 = 

n o (* + 3 ? 2 A *) - 

r = it * i, , 

g = 

3n 0 7 2 Ag , 

g" = g t + g 0 , 

h = 

3 n 0 y 2 Ah , 

h" = h t + h Q , 

K = 

n Q [l + 3r 2 ( A ^ +Ag + Ah)] , 

K" - K t + \.q . 


LONG-PERIOD TERMS 

q x = 1 - 15 9 2 q 2 = 1 - 7 6 2 
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p 3 = 

21 d 4 - 146> 2 + 1 

P 4 = 

1 

- ge 2 

P 5 = 

75 6> 4 - 30 0 2 + 11 

p 6 = 

35 - 14 6> 2 + 3 

q 7 = 

135 <y 4 - 50 e 2 + 11 

q 8 = 

45 i9 4 - 18 6> 2 + 5 

Pi = 

3e 0 2 + 4 

P 2 = 

9e 

o 2 + 4 

P 3 = 

T) 2 + Tj + 1 

P 4 = 

- 

9 17 2 - 6 17 + 7 


1 

r c J 4 

-1 



b i = ~ 4 y e o sin l o 

j 2 Pi + 5 j 2 

^2 



5 J 5 


b 2 = 7 ; q " 4 j: v pi q 3 


35 J 5 , , . 2 . 

24 To 7 e o q 4 sin x 0 


J 3 5 Js , 

J- q + 4 J- 7 P 2 q 3 


1 

4 ye 0 sin i 0 


J 4 

J 2 q s + 5 q 6 


Ja 5 Js 
j- q + 4 j: Pi L q 3 “ 6 


q 6 sm i 0 


35 Js , „ q 7 sin io 

72 To 7 6 o q 
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b n ~ T 


h 

h 


p 3 0 


T)+ 1 T 0 + 1 


1 5J5 

J q_ 4 j; 


r 2 iq 3 


Pi^ p 4 v 2 


0 + 1 T 77 + 1 


+ 9 e Q 2 + 26 


6 p 1 sin 2 i 0 
q( 0 +l) 


'10 


35 J 5 2 2 . 2 . 

72 J 2 7 e o sin x o 1^4 


3p 3 

77+I + 2 + 


_2_1 

0 + lj 


2 ^q 8 

q( 0 + 1 ) 


7 V sin i 0 

S t e - ~ (b x cos 2g" + b 2 sin g" + b 3 sin 3g" ) 


8 ii 


q 

7 e o 6 


(bj cos 2g" + b 2 sin g" + b 3 sin 3g" ) 


e 0 S 1 Z 


sin i o S x h - 


7 t) sin i ( 


(bj sin 2g" + b 4 cos g" - b 3 cos 3g" j 


ye 0 6 , 

— — (b 5 sin 2g" + b g cos g" + b 7 cos 3g" ) 


Sj \ 


7 e Q sin i 0 


(b g sin 2g" + b 9 cos g" + b 10 cos 3g" ) 


SHORT-PERIOD TERMS 


E"-e 0 sinE" = l" 


/ 77 sin 

tan t cosTT^ 


in E" \ 


e - 1 + e Q cos f" 


a - a 


0 j 1 + 2J 2 ^2 [liit 6 -7 ?) ( e2 + e7 ? + 7 7 2 ) + 3e3 sip2 ' l o cos ( 2 e" + 2f ") | 


8 2 e 


= J 2 7' 


f ffoPa 

l 1 |_ Tj + 1 


+ 3 cos f " + 3e Q cos 2 f" + e Q 2 cos 3 f " 


] 


sin 2 i o |^3 (e Q + 3 cos f " + 3e Q cos 2 f " + e Q 2 cos 3 f " ) cos ( 2g" + 2f " ) 
- t] 2 (3 cos (2g" + f ") + cos (2g" + 3f 
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J 2 y 2 0 sin i 0 |3 cos (2g" + 2f ") + e 0 [3cos (2g" + f ") + cos (2g" + 3f") 


e 0 § 2 L = - 2 J 2 y -n 


r 

-<2gj 1 (e 2 + e + t) 2 ] sin f' 


+ 3 sin 2 i , 


(? f 2 -e-e 2 )sin(2g" + f")+(e 2 + e+ 3' ! 7 2 ) sin(2g + 3f ) J 


sin i Q S 2 h = - J 2 y 2 0 sin i „ ^6 (f " - L " + e Q sin f ") 

- 3 sin (2g" + 2f") - e 0 [3 sin (2g" + f") + sin (2g" + 3f " )]| 


S 2 \ = J 2 y 2 S(l50 2 -60-3)(f"-r') + e c 


(l5 0 2 -6 0- 3) 


r k 2 + e .i 

i 77 + 1 


zzf ) 


] „ 3 n • 

v 2 - e - e 2 50 + 3 

) sin f + 2 e o s ^ n x o 

77+1 H 0 + i 


sin ( 2g" + f") 


+ 2 e o sin *0 


3 (e 2 + e + ^ t) ! 

! ) 5 0 + 3 

77+1 

+ 0 + 1_ 


sin ( 2g" + 3f " ) 


3 (50 + 3) 
+ 2 0 + 1 


sin 2 i 0 sin (2g" + 2f ")/ 


OSCULATING ELEMENTS 

The formula for the semi-major axis, a, is given above. The remaining 
osculating elements are obtained by the following: 


Be = 

8 l e + 

s 2 e , 

sin i Q 8h - sinigSjh + sin i 

Si = 

8 . i + 

S 0 i , 

e 0 8 l = e 0 8 x l + e 0 8 2 l , 
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SA - SjA + SjA, 


Mi 

- (e 0 + 8e 

) cos Z" - 

m 2 

= (e 0 + Se 

) sin l " + 

v l 

- cos h" 

sin i 0 + Si 

V 2 

= sin h" ( 

sin i 0 + Si 


e 

ii 

ti 


i 

ii 

►-■ 

o 

4- 


h 

tan“ 


l 

- t an”* 


g 

= A #/ + 


e Q S l sin l " , 

6 0 S l COS l " , 

cos i 0 ) - sin i 0 Sh sin h" , 
cos i 0 ) + sin i 0 Sh cos h" , 


+ m 2 2 • 

Si , 



Sk - h - l . 


The coordinates and velocity components are then computed from the 
osculating elements in the usual manner. 


CONCLUSION 

The formulas contained herein are valid for all eccentricities and inclina- 
tion (except inclinations near the critial inclination and i = v) yielding the same 
results (to the order of J 2 ) as Brouwer's formulas when neither the eccentricity 
nor the inclination are small. When the eccentricity is small l and g may be 
ill-defined, and when the inclination is small h, l, and g may be ill -defined. 
These cases cause no numeric problems, however, since k = g + l + h is always 
well defined. 
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